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Implications of heavy quark spin symmetry on heavy meson hadronic molecules 
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Abstract 

In recent years, many heavy mesons and charmonia were observed which do not fit in the conventional quark model expectations. Some of 
them are proposed to be hadronic molecules. Here we investigate the consequences of heavy quark spin symmetry on these heavy meson 
hadronic molecules. Heavy quark spin symmetry enables us to predict new heavy meson molecules and provides us with a method to test 
heavy meson molecule assumptions of some newly observed states. In particular, we predict an 77^/0 (980) bound state as the spin-doublet 
partner of the F(4660) proposed as a i///o(980) bound state with a mass of 4616lg MeV and the prominent decay mode i~i' c iriv. The width is 
predicted to be T(r]' c inv) = 60 ± 30 MeV. The tt + tv~ invariant mass spectrum and the line shape are calculated. We suggest to search for this 
state in B —* r)' c K 7T + 7r _ , whose branching fraction is expected to be large. 
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The chromo-magnetic interaction of a heavy quark with 
gluons is proportional to the magnetic moment of the heavy 
quark, and is suppressed by the heavy quark mass (for re- 
views of heavy quark symmetry, see Refs. HI 0])- There- 
fore in the heavy quark limit mq — > 00, the interaction is 
spin-independent, and a new symmetry appears called heavy 
quark spin symmetry. Due to this symmetry, there are spin 
multiplets of both heavy mesons and heavy quarkonia, as e.g. 
the {D,D*} and {rj c ,J/ip}, The masses of the members 
within the same spin multiplet would be degenerate in the 
heavy quark limit. In this Letter, we extend this symmetry 
to possible heavy meson molecules, which were observed in 
recent years (for reviews, see e.g. Ref. JslHIslEl)- By heavy 
meson molecules, we mean bound states consisting a heavy 
meson/heavy quarkonium and a light hadron, or two heavy 
mesons. In this Letter, we will focus on the former type. 

As an example, let us focus on the D* (2317) and the 
r> s i(2460) first. They were proposed to be S-wave hadronic 
molecules whose components are mainly DK and D*K, re- 
spectively J3, HI IU QH U , 12 ] (for the latest development on 



the J * n (2 317) in the hadronic molecular picture, we refer to 
Refs. jrUm HO])- Their masses are measured to be |QT 



M 
M 



D* (2317) 
D s i(2460) 



2317.8 ±0.6 MeV, 
2459.6 ±0.6 MeV. 



(1) 



Were they the bound states of DK and D*K, respectively, the 
binding energies are 



e D* (2317) — M]j 
£D s i(2460) = ^Id' 



H M K - M D » o ( 23 i7) = 
+ Mk - -Md s1 (2460) 



45 MeV, 
= 45 MeV, (2) 
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where we have taken the averaged masses within the same 
isospin multiplets of D, D* and K. One notices that the 
binding energies are the same. For molecular states this ap- 
pears to be natural: first of all the leading interactions of 
light mesons with D and D* mesons are independent of the 
heavy quark spin and secondly the light meson-Z^*' me- 
son Greens functions, which provide an important input to 
the bound-state equations, are to a very good approximation 
mass-independent as long as evaluated close to the corre- 
sponding threshold. As a result, hadronic molecules also fall 
in spin multiplets, and, most importantly here, the splitting 
within one multiplet remains the same as the hyperfine split- 
ting between the heavy mesons which are the components of 
the hadronic molecule. 

In the same way the hyperfine splitting within a heavy 
quarkonium spin multiplet will also be untouched by the in- 
teractions with light mesons. The interaction between a heavy 
quarkonium and light hadrons occur mainly through exchang- 
ing soft gluons. The leading order heavy quarkonium inter- 
action with a soft gluon field comes from the chromo-electric 
dipole interaction 1 18 , lH 2(J which is spin-independent. The 
chromo-magnetic interaction is suppressed by I/tuq l2~lll . 
Since any hadron should be color singlet, the number of the 
exchanged soft gluons should be at least two. Therefore, 
the suppression of the spin-dependent interactions between a 
heavy quarkonium and light hadrons is at least 1/toq. As a 
result, a bound-state of a heavy quarkonium and light hadrons, 
called hadro-charmonium in Ref. 12211 . will have partner(s) 
whose components are the same light hadrons and the spin- 
multiplet partner(s) of the same heavy quarkonium. The mass 
splitting within the molecular spin-multiplet will be, to a very 
good approximation, the same as the heavy quarkonium hy- 
perfine splitting. 

This nice feature enables us to predict new heavy meson 
molecules and provides us with a method to test heavy me- 
son molecule assumptions of some newly observed states as 



1 



illustrated in the following. 

The Belle Collaboration [23] observed a resonant structure, 
called F(4660), in the i[;'ir + ir~ final state using the method 
of initial state radiation. The line shape of the state was fitted 
with a P-wave Breit-Wigner 13011 in the experimental paper, 
and they got 4664 ±12 MeV. In Ref. we demonstrated 
that the experimental data support a tp'f q(980) bound state 
hypothesis for the y(4660). As a result of fitting the mass, 
which gave M Y = 4665±f MeV, we could calculate the spec- 
tral distribution. If this interpretation of the Y (4660) is indeed 
correct, heavy quark spin symmetry implies that there is an 
??c/o(980) bound state, to be called Y n in the following. The 
quantum numbers of such a state are J p = — . The <q' c and 
ip' lie in the same spin multiplet, and their mass splitting is 
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FIG. 1 
decay. 
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The 7r + 7r _ invariant mass spectrum for the Y n 



i + - 



AM = M$i - M. 



49 ± 4 MeV. 



(3) 



From the above analysis, the mass of the rj' c fo (980) bound 
state would be 



My 



M 



Y"(4660) 



AM = 4616j1 MeV , 



(4) 



where the uncertainties of the mass of the F(4660) and AM 
were added in quadrature. Similar to the Y(4660) decaying 
predominantly into ip'-mr, the dominant decay channel of the 



Y v would be Y v 



r] c Trn. 



If two particles form an S-wave bound state which is very 
close to the threshold, there is a way to model-independently 
connect the effective coupling constant of the bound state to 
its constituents, g, directly to the molecular admixture of the 
state 11251 12611 . Historically, Weinberg used this method to 
show that the deuteron is not an elementary particle. Espe- 
cially, one may write for a pure molecule 



9 2 A , ^2 /2e 

— = 4(mi + m 2 ) W — 

47T V M 



1 + 



/3 



(5) 



where mi and m-i denote the masses of the constituents, e the 
binding energy related to M, the mass of the molecule, via 
M = mi + m-2 — e,/i = mim2/(mi + 7712) the reduced 
mass, and 1/(3 the range of the forces. 

With the effective coupling constant fixed by Eq. (0, we 
can predict the irir invariant mass spectrum and the decay 

Denoting the tt + tt~ invariant 



width of the Y, 



n 



r]' c ir + ir 



mass by m n „, we have for the differential width 



g 2 q 



8ttM2 

1 7 



[7T + 7T 1/ \ 

■P fo 



(6) 



where q is the magnitude of the three-momentum of the r]' c in 
the Y n rest frame 



My 



(m^ + M V ,J' 



M Y 

1 11 



(m nv - M^y 



and 



2Mv 



Pf (totttt) is the 7r + 7r fraction of the /o spectral 



function, 



(m n7V ) = - 

7T 



.in 



(7) 



where ft^(m w ) = U.f p (m^) ~ Re(nj o b (m /o )) denote the 
renormalized self-energies of the /o with respect to the chan- 
nel ab = 7T7T or KK. Analytic expressions are given in 
Ref. 12711 . The imaginary part of the self-energy of the /o 



is fixed by unitarity 

Im ( n /o +,r ( m ^)) = m /o r /o->?r+^-( m ™)' 



and 



4m£ 



167T771 



(8) 



(9) 



fa 



The input parameters related to the /o(980) are taken 
from the fits provided in Ref. [28 J. To be specific, we use 
m/ = 0.9862 GeV, gf K+K- = 3.87 GeV, and g :fo1T + 7T - = 
—2.03 GeV, which are the central values of the various param- 
eters of fit K2 shown in Table 4 of that reference. The cou- 
plings for the neutral channels are fixed using the isospin rela- 
tions. With these parameters, and M^ = 3637 ± 4 MeV lfl7ll . 
the 7r + 7r~ invariant mass spectrum for the Y„ decaying into 
i]' c ir + ir~ is predicted in Fig.Q] where the solid line shows the 
result using the central values of M J( < and My v ■ 

Integrating over the irir invariant mass, and considering 
both the r)' c ir + ir~ and rj' c ir°ir° channels, we get a width of 
the Y v — > r]' c Trir as 58 ± 5 MeV. Here we only took the lead- 
ing term of the effective coupling constant as given in Eq. ([5]), 
and the uncertainty comes solely from that of the masses of 
the Y v and r/' c . However, to get a more realistic estimate of the 
uncertainty, we need to also estimate, e.g., higher order terms 
of the effective coupling constant. As stated before, the in- 
teractions between a charmonium and light hadrons are medi- 
ated by soft gluons, so we may estimate the range of forces as 
1/(3 sa 1/Aqcd- Noticing y/ifle « 100 MeV, the uncertainty 
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FIG. 2: Line shape of the Y v in the t]' c -k + tv~ invariant mass distribu- 
tion. For clarity we only show the line corresponding to the central 
values of the parameter space. Note that the asymmetry seen in the 
distribution is an unavoidable consequence of the molecular struc- 
ture. 



Py c7! *" ' (M). That is the line shape of the Y n in the rj' c ir + ir 
mass distribution as given in Fig. [2] 

The proposed ?y^/o(980) bound state can be searched for in 



B decays. We suggest to search it in B d 



Taking data from Ref. lfl7ll for three measured channels, the 
branching fraction of the so far unmeasured decay B ± — > 
rj' c K ± 7r + 7r~ can be estimated as 



B(B ± -» 
= B(B ± 



r]' c K ± TT + n 



(3.4 ± 1.8) x 10" 



1 x 10" 



B(B ± -> iP'K ± Tr+ir-) 

B(B ± -> ijj'K*) 
i (1.9 ± 1.2) x 10" 3 
(6.48 ±0.35) x 10" 4 



(13) 



of g from this source, and hence the width, is about 50%. 
We want to stress that this estimate is clearly quite conserva- 
tive. Taking the expected inverse mass of the lightest glueball 
as the range of forces might appear equally justified (and this 
would lead to a much smaller uncertainty). On the other hand 
we do not explicitly include uncertainties from other sources 
like higher orders in the 1 jniQ expansion. Thus, as result, we 
get 



r(Y n 7/>tt) = 60 ± 30 MeV 



(10) 



and the uncertainty in the ir + ir~ invariant mass spectrum is 
reflected as the band in Fig.Q] The uncertainty in the coupling 
constant has a much larger effect on the signal than the one of 
the mass of the Y n . 

Assuming the width of the Y n being saturated by the r)' c mr 
final state, the line shape of the Y v can also be predicted. For 
this we use a dispersion integral, which gives us an expression 
for the Y self-energy, Ily (M), for arbitrary values of M 



n y „(M) = i 

7T 



ds — 

M 2 s 



MyT^Hyfs) 



AP 



ze 



(11) 



where M t h r = M v ' +2m w denotes the lowest physical thresh- 
old of relevance here and Ty* (y/s) the total width of the Y v as 
a function of y/s. Note that this treatment is completely con- 
sistent to what was done for the fa. With the self-energies at 
hand we may now give the expression for the spectral function 
of the y„ (4616) 



MyT^HM) 



M 2 



My 
XT) 



(12) 



where, as above, we defined fly (M) = Ily (M) — 
Re(ITy i) (My i) )). Replacing the total width of the Y v in 
the numerator of Eq. ( fl~2l ) by T(Y ri — > 7f c ir + ir~), one gets 



Such a large branching fraction offers a great opportunity of 
finding the Y v in the B decays, although we cannot predict 

B{B ± -> Y n K ± ). 

In summary, the heavy quark spin symmetry, which is ex- 
act in the heavy quark limit, is extended to the systems made 
of a heavy meson/quarkonium and light hadrons. We argue 
that the hyperfine splitting remains untouched in heavy me- 
son molecules. Based on this observation, there should be an 
7^/o(980) bound state with a mass of 4616^g MeV, were the 
y(4660) a V7o(980) bound state as suggested in Ref. 0]. 
Such a bound state would decay mainly into i]' c inr with a 
width of 60 ± 30 MeV. In addition, analogous to the Y (4660), 
we also predict decays into rj' c K + K~ and 77^77. There is also 
the possibility of a decay into A+A^. We predict the 
invariant mass spectrum of the Y^(4616) decay into rj' c 7r + Tr~, 
and the prediction is parameter-free. We also predict the line 
shape of the state in the t]' c / k +/ k^ final state assuming its de- 
cays are saturated by the rj' c irir. 

The state can be searched for in the S-factories Belle and 



BaBar. The branching fraction of the B 



± 



7T 7T IS 

estimated to be of order 1 x 10~ 3 , and hence there is a great 
opportunity to find the 77^/0 (980) bound state proposed here 
in the rj' c ir + ir~ final state of this d ecay . The state can also 
be studied with PANDA at FAIR Q in the future. Such 
a study would be helpful to understand better not only the 
XYZ states observed in recent years but also the interaction 
between a charmonium and light hadrons, which can provide 
useful information for understanding the charm production in 
relativistic heavy ion collisions. 
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